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Wave propagation in nonlinear theories of the electromagnetism described by Lagrangian densities 
dependent upon its two local invariants L{F, G) is revisited. On the light of the recent findings in 
metamaterials, it is here shown that trirefringence is also a possible phenomenon to occur in the 
realm of such nonlinear theories. A specific model exhibiting this effect is investigated both in terms 
of phase and group velocities. It is claimed that wave propagation in some well known nonlinear 
models for spin-one fields, like QED and QCD in certain regimes, may exhibit trirefringence. 
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I. INTRODUCTION 

As it is well known, nonlinear theories of electromag¬ 
netism exhibit birefringence phenomenon. The most 
popular example appears in the quantum electrodynam¬ 
ics (QED) where polarization effects are activated in the 
limit of large fields ~ E,,,. = /eh), inducing an 

effective optical axis in the vacuum. In such situation, 
a light ray is expected to ^lit in two rays propagating 
with different velocities [ll-lj . An experimental setup de¬ 
signed to measure the birefringent properties of the QED 
vacuum was long ago proposed However, direct mea¬ 
surements of this effect are not yet conclusive and are still 
under consideration Q. The influence of the nontrivial 
vacua on the propagation of electromagnetic waves was 
discussed in several distinct physical configurations 0- 
fl^ . Conditions for the occurrence of birefringence of 
gluon fields was also studied [T^ . 

In the context of material media, birefringence effects 
are expect to occur in several distinct situations. It oc¬ 
curs naturally, for instance in certain crystals present¬ 
ing optical axes min, or artificially when optical axes 
are induced b y m eans of external applied electromag¬ 
netic fields [ia,[i3- Nowadays, birefringent materials and 
methods including this effect have been incorporated in 
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several technological devices M- Birefringence is also 
a powerful optical tool to investigate properties of new 
materials, biological systems and others [l9l - l2ll| . 

On the other hand, trirefringence was only recently 
considered as a possible phenomenon in material media. 
It was measured ^ in tailored photonic crystals [^ . 
and the theoretical description of this effect in media 
characterized by effective dielectric coefficients was pro¬ 
posed [2^. In this case [2^, only when some of the di¬ 
electric coefficients are ne gat ive, could trirefringence take 
place. Metamaterials [251 - 12^ seem to be good candidates 
for supporting this effect, due to the controllability of 
their dielectric tensors. With the present day technology 
of producing such new media, it is expected that trire¬ 
fringence will play some important role in technology of 
optical systems, as birefringence has done. 

Usually, effects occurring in the realm of Maxwell elec¬ 
tromagnetism in material media are expected to occur 
in the realm of nonlinear electromagnetic theories. It 
is possible to build up analogue models between these 
two domains where the coefficients describing a specific 
dielectric medium are mapped as derivatives of the La¬ 
grangian density describing a nonlinear theory. In this 
way, trirefringence should also be a possible effect in non¬ 
linear electromagnetism. By deriving and using the gen¬ 
eral description for wave propagation in the limit of geo¬ 
metrical optics, in this paper we show that trirefringence 
is in fact a possible effect in nonlinear electrodynamics. 
It is not our purpose to set the general conditions a model 
must fulfill in order to present trirefringence, but only to 
show the effect as a possible one in the domain of nonlin- 
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ear electromagnetism. A particular model is thus inves¬ 
tigated where such phenomenon is shown to occur pro¬ 
vided that convenient external fields are set. The phase 
and group velocities of the waves are derived, as well as 
the corresponding polarization vectors. A numerical ex¬ 
ample is graphically studied. 

The model examined in the paper corresponds to the 
effective Lagrangian density for QED in the regime of 
large fields. The trirefringence phenomenon is shown to 
occur whenever the model applies, although its measura¬ 
bility requires the control of very large fields. A possible 
arena to search for this effect could be the special fluids 
recently produced by high energy collisions, as addressed 
later in the concluding section. Further, this model is 
also useful in the context of analogue models in mate¬ 
rial media, where the predicted effect could be tested on 
optical systems, for instance in metamaterials. 

In the next section, nonlinear electromagnetism is 
briefly revisited and the field equations are presented in 
terms of the general two parameters Lagrangian density. 
In Sec. uni the corresponding wave propagation is ex¬ 
amined. The eigenvalue problem is stated and solved, 
resulting in the general fourth degree equation for the 
phase velocities. This equation is solved for a specific 
nonlinear model in Sec. El The corresponding polariza¬ 
tion states and the description of the effect in terms of 
group velocities are also discussed. Conclusions and final 
remarks are presented in Sec. 0 

Throughout this paper we employ the Minkowski met¬ 
ric = diag(-|-I, — I, — I, — I). The completely skew- 
symmetric tensor is defined by = I. We set 

the units such that the velocity of light in empty space 
is c = 1. 


II. NONLINEAR ELECTRODYNAMICS: FIELD 
EQUATIONS 

Nonlinear Abelian theories for electromagnetism can 
be formulated by means of the general Lagrangian den¬ 
sity L = L{F,G), where F and G are the two local gauge 
invariants of the electromagnetic field. These invariants 
are defined in terms of the electromagnetic tensor field 
Ffj_,y, and its dual 

FaP = -jVap''''Far, ( 1 ) 

as 

F = Ff^'^Ff,, (2) 

G = F^".F^,. (3) 

In terms of the electric E and magnetic B field strengths 
we have F = —2{E^ — B^) and G = —AE ■ B. 

The field equation can be obtained from the least ac¬ 
tion principle and it can be presented as Q, 

= 0, (4) 


where is defined by 

^ LFFF^'^F°‘f^ + LggF^''F°^^ 

+Lfg (^F>^''F°‘^ + Ff^''F°‘l^'^ . (5) 

Use is being made here of the notation = 

d'^L/dX^dX'^ ■ ■ ■ dX"^ previously introduced Q, where 
each X® is one of the two invariants E or G upon which 
the Lagrangian L arbitrarily depends. We notice that 
the above defined rank-4 tensor presents the following 
symmetries: = -iV'®®''^" and 

In addition to Eq. o, F^i, satisfies the Bianchi idem 
* 

tity E^'^= 0, which implies in the existence of a poten¬ 
tial vector as 

F^u = - Avf- ( 6 ) 


III. NONLINEAR ELECTRODYNAMICS: WAVE 
PROPAGATION 

Let us now discuss the propagation of electromagnetic 
waves in the general formulation of nonlinear electro¬ 
dynamics. We restrict ourselves to the propagation of 
monochromatic waves in the limit imposed by geomet¬ 
rical optics Bin. The method of field discontinuities 
will be used, which can be briefly stated as follows [^1^. 

Consider a differentiable inextendible oriented border¬ 
less hypersurface E, defined locally by = 0, where 

(/) is a real differentiable scalar field which locally is a 
function of the spacetime coordinates = {t,x). Let 
be the spacetime points whose coordinates satisfy 
(j){x^) > 0, and similarly U~ be such that (j){x^) < 0. 
Let P be any given point of E. For each sufficiently 
small r > 0, let Vr{P) be a neighborhood of P which 
consists of the spacetime points Q whose Euclidean dis¬ 
tance from P is [(fg — -|- ||xq — smaller 

than r. Let P+ e [/+ n Vr{P) and P~ G U~ fl Vr{P) be 
any two neighbor points from P arbitrarily chosen at op¬ 
posite sides of E. Let / be any given tensor field defined 
at Vr{P). The Hadamard discontinuity at P of / across 
E is defined as 

[/]s(P)= lim [/(P+)-/(p-)]. (7) 

r—>^0+ 

Suppose / such that [/]e = 0 for each P e E. Follow¬ 
ing Hadamard [l^, we have [/,a]s(T’) = k\f{P), where 
k\ = 4',\\p is the normal vector to E at P and / is a ten¬ 
sor field defined at E with the same rank and the same 
algebraic symmetries as those of /. 

We assume the electromagnetic tensor field P^^ to be 
smooth in each P®*®, but merely continuous at E (that 
is to say, the F^^ are continuous functions at E but 
their derivatives may present discontinuities at E). The 
Hadamard discontinuities at E of Eq. and of the 
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derivative of Eq. dd]) lead to Q 

fpxk^ + = 0 (8) 

Lf 

and 

fajH — ^akj^ Ci^kiy, (9) 

where the quantities fa/s are related to the derivatives of 
Fai 3 on E by [Fap^xlx: = fapk\ and is the polariza¬ 
tion vector [Ap^ap]^ = Sfikakis. We set k\ = ujV\ + q\ 
as the wave 4-vector, where V\ = is the 4-velocity 
of the observer which decomposes F^^'' into electric and 
magnetic fields. The components of this 4-vector kx are 
thus the frequency w and the wave vector q = qq, where 

ll<?P = = 1- 

Taking together E™. JS]) and ([5]) , we obtain the general 
eigenvalue equation [ll, Q 

= 0 , ( 10 ) 

where 

( 11 ) 

Lp 

with = k^k\. Nontrivial solutions of Eq. (US can be 
found only if det | Z^i, | = 0, the well known generalized 
Eresnel equation, and yields 

= ( 12 ) 

where p = F°‘>^Fa'^kpkv, and 

a = Pp + ‘iLpiCLpc — FLcg) 

-{LffLgg - (13) 

P = ALf{Lff F Lgg) — 8{LffLgg — L%g)^^ ( 14 ) 

'y = 16 {LffLgg - L%g)- ( 15 ) 

One can also recast the quantity p that appears in Eq. 

(HH) as 

p = Ef -PE^ + {q-Bf -q^B^ + 2ujpE X B. (16) 

The phase velocity v = oj/q of the electromagnetic 
waves can be obtained from Eq. (HD). In fact, it is 
straightforward to show that this equation can be pre¬ 
sented as a fourth-degree equation for the phase velocity 


V as 

041 '"' -I- 03 ?;^ -|- a2V^ + aiv -|- oq = 0, (17) 

where we have defined 

a4, = a - PE^ + jE^ (18) 

03 = 2(/3 - 2-fE^)q -ExB (19) 

02 = -2a + p[E^ -B^ + {P Ef + {q ■ Bf] 
+2-f{2{q -ExBf- [(g • Ef + {q ■ Bf 
-B^]Ep ( 20 ) 

a, = -2{p-2Piq-Ef + {pBf 

-B^]}PExB (21) 

oo = a + p[B^ - (q ■ Ef - [q ■ Bf] 

F^[B^-{q.Ef-{q-Bff. (22) 


As stated by Eq. (El), we can find up to four solutions 
for the phase velocity in the same wave direction. In 
the next section we will analyze some special cases where 
multirefringence phenomena may occur. 

Dispersion relations for light propagation in nonlinear 
electrodynamics can also be investigated by means of the 
photon mass operator Bmi. In such context the prop¬ 
agation of photons in homogeneous magnetic field was 
investigated long ago [s^ and birefringence phenomena 
was described for some field configurations. 


IV. A MODEL FOR TRIREFRINGENCE 

In what follows we shall study a particular model for 
nonlinear electromagnetism which presents interesting 
multirefringence features. Let the one-parameter non¬ 
linear Lagrangian density 

Lnl =-^hoFXog^, (23) 

where &o and A are constants. Particularly these con¬ 
stants can be chosen in order to split the above model in 
a Maxwellian part plus a nonlinear contribution. 

This model appears in different contexts in the liter¬ 
ature, for instance as the effective Lagrangian density 
for quantum electrodynamics (QED) [33| in the regime 
of large fields . In this case the constants in Eq. 

dMl) are related to the fine structure constant a and to 
the critical electromagnetic fields for which vacuum po¬ 
larization effects begin to become important. Further¬ 
more, with an appropriate choice of dielectric coefficients 
this model can describe several kinds of magnetic materi¬ 
als which can be used as analogue systems to investigate 
properties of the vacuum of the non-Abelian gauge field 
[38]. A more detailed discussion about this issue is pre¬ 
sented in the concluding section. 


A. Phase velocity 


The study of wave propagation in this special case can 
be done by following the lines presented in Sec. IIIII We 
will investigate multirefringence phenomena using Ljyl 
in the Abelian case only. 

To study a simplified case, let us assume constant ex¬ 
ternal electric E = Ex and magnetic B = By fields, 
much larger than their wave counterparts. The wave vec¬ 
tor q is assumed to lie in the xz-plane so that, p y = 0, 
q ■ X = sin0, and p z = cos 9. Thus, 6 is the angle 
between q and z directions. From the above notation, 
p E X B = EB cos9. Taking Ljvl in Eq. (flTll we obtain 
the following results 


^^0 = 1, (24) 

2EBcos9± x/(x-2B2)(x-2EPos^9) , , 

v± =-,(25) 


X 




where it was introduced the shortcut 


4 


X = 2E^ -^(l + log . (26) 

The quantity vg is isotropic and does not depend on any 
choice for the configuration of the fields nor direction of 
propagation. However, the other solutions will be dif¬ 
ferent for different configurations of fields and q direction, 
set by the angle 0. As one can see, three different veloci¬ 
ties in the same direction occur provided the square-root 
is smaller than 2BEcosO. This naturally imposes con¬ 
ditions on the fields. We observe that trirefringence will 
occur in a region defined hy —9c < 0 < 0c^ where 

0c = arccos | (s + log ^ | . (27) 

Therefore, it does exist iff 

Q<F <Fc, = e-3 ~ 0.0498 (28) 

and 

E>Em, A^ = -|(3 + log J). (29) 

Birefringence takes place in the regions 0c < 0 < — 0c 

and —(tt — 0c) < 9 < —9c- For the remaining angles, just 
the ordinary wave exists. When F = Fc, trirefringence 
is not present for any direction. If one keeps F fixed and 
satisfying Eq. (1^ and decreases the electric field [satis¬ 
fying Eq. ([29|) ] , then the region where trirefringence takes 
place decreases and the difference between the moduli of 
the extraordinary solutions increases. In Fig. [I] the nor¬ 
mal surfaces associated with Eqs. (|24ll and (|25ll 

are plotted for a selected set of the fields. It can eas¬ 
ily be seen there the region exhibiting trirefringence, as 
anticipated, for fields satisfying Eqs. (l28l) and (l29l) . 


B. Polarization 


Let us now return to the Fresnel-like eigenvalue Eq. 
OT. The matrix given by Eq. m reduces, for the 
Lagrangian in Eq. ([23|) , to 



- m + 

+ iq X B)^ 


vEf^ -h (q ■ E)V>^ 


F(l + log^) L 
vE^ -I- {q ■ E)V„ + {q X B)i, . (30) 


This suggests that the polarization vector e’' should con¬ 
veniently be decomposed as a linear combination of the 
three vectors which appear in as 

= avE'' + b{q ■ E)V'' + c{q x B)'' + dk‘', (31) 


Vx 



FIG. 1. (color online). Normal surfaces for the nonlinear theory 
described by Eq. (I23II . The ordinary wave is represented by the 
circular thick line and the extraordinary waves are represented by 
the dashed and dot-dashed curves. The symbols -I-, and — and 
o indicate the solutions presented in Eqs. (12511 and (12411 . The 
chosen values of the electric and magnetic fields are B = O.IA 
and E = 0.09A. These values satisfy Eqs. (12811 and (12911 : hence 
trirefringence is present. Since this is the case, also birefringence 
and one refraction must take place. The regions where these 
effects take place are related to the critical angle, 9c, as given 
by Eq. (I27II . The trirefringent region lies between the dotted 
straight lines, while the birefringent regions are constituted by 
the angles limited by the dotted and dot-dashed straight lines. 
The region presenting one refraction lies between the dot-dashed 
straight lines. 


where a, b, c are arbitrary constants with the same phys¬ 
ical dimension. The fourth term, with an arbitrary con¬ 
stant d, was introduced because Eq. m remains un¬ 
changed by it. Equation OT then reads 

[a(t>2 - 1) - Y]vE>^ + [biv^ -1)-Y]{q- E)V>^ 

+ [c{v^ - 1) - r] (9 X H)^ = 0, (32) 


where E is a shortcut for 


Y = 


4 

-|-c(9 X B)^ - 


av^E^ - 
b{q-E)^' . 


(a -I- c)vq ■ {E X B) 


(33) 


For the = 1 case, Eq. ([^ yields E = 0, from 
which the polarization state is given by a = {q ■ B)^, b = 
E^ — {qx B)^, c = -{q- B)^ up to a global multiplicative 
factor. For the yf 1 case, Eq. (15^ yields a = b = c = I 
up to a global multiplicative factor. Assuming E ■ B = 0 
and q ■ B = 0, then the phase velocities for this case 
are given by Eq. dUD. Once the physical parameters 
a, b, c were found in either case, then the particular gauge 
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choice d = —b{q ■ E)/{qv) ensures e'' to lie in the space 
orthogonal to V. 


C. Group velocity 

It is well known that, in the geometrical optics approx¬ 
imation the wave equation is a linear equation for the 
perturbed fields, even in the context of nonlinear electro¬ 
dynamics [ij. In this case wave packets can be build up 
by superposing plane wave solutions, whose phase veloc¬ 
ities were obtained above. Thus, it is important to deal 
with the group velocities u = Aw/Aq for the propagation 
analysis [iJ. For the particular case of the nonlinear 
Lagrangian density iwL, with the same configuration of 
fields and wave vectors we assumed above, the associated 
group velocities are 


U = UxX -\- UzZ, 


(34) 


where Ux = sin 9 and Uz = cos 9 for the ordinary u = 1 
mode, thus stating that the group and phase velocities 
coincide for this case. For the extraordinary modes, we 
obtain 


I yui — AEBv± cos 9 -I- 2E‘‘ cos"^ 9 
= - „ „„ --- sin0, 


Xv± — 2EB cos 9 


(35) 


and 


X'*^± cos 9 — 2EBv± cos 29 — 2E^ cos 9 sin^ 9 
^ — 2EB cos 9 


(36) 

where we are considering v± as given by Eq. (1251) . Hence, 
two extraordinary rays and one ordinary ray can be found 
in the xz- plane. But their dependence on direction must 
be investigated more carefully, as we shall do in what 
follows. If one defines (p as the angle between the group 
velocity and the z-axis, then (p = 9 ior the u = 1 case. For 
the extraordinary modes, it follows from Eqs. (ESI-dMl) 


that 


(X'c± “ AEBv± cos 9 + 2E^ cos^ 9) sin 9 

tan p = — ---=—, 

Xu± cos 9 — 2EBv± cos 29 — 2E^ cos 9 sin^ 9 

(37) 

which gives us as a function of 9. For the case of the 
extraordinary solutions, the analytical inversion of Eq. 
(1371) to give 0 as a function of p is very involved. Hence, 
numerical analyses turns out to be more clarifying. Eig. 
[2] summarizes such a numerical analysis for the same pa¬ 
rameters assumed in Eig. [T] One notices that the region 
of the xz plane where the ordinary and extraordinary 
group velocities can be found constitutes a trirefringent 
region. For the complementary region of the plane of 
propagation, there exists just the ordinary group veloc¬ 
ity. 

Comparing Figs. [T] and [2] we find angular sectors for 
which plane waves associated with the extraordinary po¬ 
larization modes are supposed to propagate, but with 


Ux 



FIG. 2. (color online). Ray velocities for the nonlinear model 
described by Eq. (I23II . The fields are the same as in Fig. [T] 
The symbols -I-, — and o indicate the solutions presented in Eqs. 
(12511 and II24II when substituted in Eqs. II34I) - (I36I) . As it can be 
seen, trirefringence occurs in the region lying between the dotted 
straight lines. In the complementary region only the isotropic ray 
solution propagates. 


no propagation of the corresponding wave packets. This 
feature relies on the fact that, in our analysis, the wave 
packets for the extraordinary modes do propagate along 
directions which are not generally equal to the directions 
of the corresponding plane waves components, as it is ex¬ 
plicitly shown in Eq. dSZl)- This is due to the fact that 
the phase velocity is dependent upon the direction of the 
wave vector; hence w = v{q)q. This dependence leads to 
a term in the group velocity that is perpendicular to the 
phase velocity. The magnitude of such a term may be 
comparable to the magnitude of the phase velocity itself, 
yielding therefore to a possibly different behavior of the 
two aforementioned velocities. 


V. CONCLUSION AND DISCUSSION 

Trirefringence phenomena is not an effect exclusively 
occurring in nonlinear metamaterials [2^ [2^ . As shown 
here, it is possible to formulate a nonlinear model describ¬ 
ing electromagnetism where this effect is also expected to 
occur. Possible extensions of the presented model can be 
sought by adding the dual invariant G, or else trying 
other nonlinear Lagrangian models. 

As it is well known, in the regime of small fields QED 
is governed by the Euler-Heisenberg effective Lagrangian 
density [s^- In this situation birefringence effect is pre¬ 
dicted to occur [l| . Experiments are still under considera¬ 
tion in order to confirm such prediction im. However, 
when the regime of large fields is considered, the effec- 











6 


tive Lagrangian governing QED presents the same form 
as the nonlinear model stated in Eq. (IMl), leading to the 
conclusion that trirefringence phenomenon is expected to 
occur in this regime. 

Multirefringence phenomena could also be found for 
systems described by nonlinear Lagrangian densities 
which depend on non-Abelian gauge fields. In such cases, 
the field strength tensors would not be gauge-invariant. 
Nevertheless, it can be easily shown that the propaga¬ 
tion of the field disturbances would be described by the 
same equations presented above, since only second-order 
derivatives of the gauge vector field may present non-zero 
Hadamard discontinuities. 

The nonlinear model discussed in the text was pro¬ 
posed long ago [MS mi as the effective Lagrangian 
of quantum-chromodynamics (QCD) or other Yang-Mills 
theories with non-trivial vacuum properties. When con¬ 
sidered for this purpose, the Lorentz invariant parameter 
F in L{F) is extrapolated to be where the 

index (a) runs in the inner Non-Abelian group space. 
Considering the regime of small coupling and taking the 
limit of large mean fields {F/\^ ^ 1) a Lagrangian den¬ 
sity with the same functional form as appearing in Eq. 
(| 2 ^ is obtained [M] • In this context 69 may be identified 
as a /3-function coefficient at leading order and A a con¬ 
stant related to the mass scale. It was claimed [M| that 


the form for Latl may also result in the leading terms 
in the limit of weak fields {F/X^ 1), as in both cases 

|log(F/A^)| 1. Hence trirefringence of non-Abelian 

gauge fields can also occur and in principle it can be 
observed, provided the required field configuration is ap¬ 
proached. The practical arena for such kind of study is 
a quark-gluon plasma (QGP), recently observed in high- 
energy heavy ion collision experiments where the gluon 
field is deconfined and can propagate in the bulk of the 
QGP. Some symmetric field configuration has been in¬ 
vestigated, and a possible observable has been proposed 
[ 1 ^ . With the progress of measurements on asymmetries 
in the experiments, which is now a hot topic in RHIC 
and LHC, other kind of field configurations can be fur¬ 
ther investigated. 
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